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Abstract. We investigate the properties of the standard perturbative expansions which describe the early stages 
of the dynamics of gravitational clustering. We show that for hierarchical scenarios with no small-scale cutoff 
perturbation theory always breaks down beyond a finite order q+ . Besides, the degree of divergence increases with 
the order of the perturbative terms so that renormaiization procedures cannot be applied. Nevertheless, we explain 
that despite the divergence of these subleading terms the results of perturbation theory are correct at leading 
order because they can be recovered through a steepest-descent method which does not use such perturbative 
expansions. Finally, we investigate the simpler cases of the Zel'dovich and Burgers dynamics. In particular, we 
show that the standard Burgers equation exhibits similar properties. This analogy suggests that the results of the 
standard perturbative expansions are valid up to the order q+ (i.e. until they are finite). Moreover, the first "non- 
regular" term of a large-scale expansion of the two-point correlation function should be of the form R~ 2 a 2 (R). 
At higher orders the large-scale expansion should no longer be over powers of a 2 but over a different combination 
of powers of 1/R. However, its calculation requires new non-perturbative methods. 
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1. Introduction 

In standard cosmological scenarios large-scale structures 
in the universe have formed through the growth of small 
initial density perturbations by gravitational instability, 
Besides, 



see 



Peebles (1980) 



in most cases of cosmologi- 
cal interest the power increases at small scales as in the 
CDM model ( |Peebles (1 982)1 ) . This leads to a hierarchi- 



cal scenario of structure formation where smaller scales 
become non-linear first. Therefore, one often uses a per- 
turbative approach in order to describe the density field 
at large scales or early times. This is of practical inter- 
est in order to study for instance weak lensing effects 
on large scales, which may in turn constrain the cosmo- 
logical parameters (e.g., Bernardeau et al. (1997)). The 
standard perturbative expansions are derived from the 
equations of hydrodynamics, coupled with Poisson equa- 
tion, by writing the density field as a perturbative ex- 
pansion over the initial density contrast (or rather over 
the linear growing mode). Such a procedure is detailed 
in [Fry (1984)| and poroff et al. (1986)1 As noticed in a 



previous work (paper 1) one obtains the same perturba- 
tive expansion from the collisionlcss Boltzmann equation. 
This clearly shows that perturbation theory cannot handle 
shell-crossing (in fact, it does not capture any effect aris- 
ing from multi-streaming) and it only provides asymptotic 
expansions. Moreover, as noticed in Lokas et al. (1996) or 



[Scoccimarro fc Frieman (1996) the terms encountered in 
perturbation theory actually diverge at next-to-leading 
order for steep linear power-spectra -Pl(^) oc k n with 
n > — 1. This shows that perturbation theory is rather 
ill-defined and it casts some doubts on the validity of its 
predictions, even when they are finite. 

Thus, in this article we reconsider the properties of the 
standard perturbative expansions which describe the dy- 
namics of gravitational clustering. First, we recall in Sect J| 
the derivation of these perturbative expansions from the 
equations of motion. Then, we show in Sect.||that for hi- 
erarchical scenarios where the power increases at small 
scales with no bound (e.g., linear power-spectra with n > 
—3) perturbation theory always breaks down beyond a fi- 
nite order q + as divergent terms appear. However, we ex- 
plain in Sect.|| that, despite these divergences, the results 
of perturbation theory for the cumulants of the density 
contrast are correct at leading-order because they can also 
be obtained from a steepest-descent method which does 
not make use of such perturbative expansions. Finally, we 
consider in Sect .[3] the simpler cases of the Zel'dovich and 
Burgers dynamics. In particular, we show that the stan- 
dard Burgers equation shares many properties with the 
gravitational dynamics. This suggests that the results of 
perturbation theory are valid up to the order q + (i.e. until 
they are finite) and that beyond this order the large-scale 
expansion of the two-point correlation (for instance) is no 
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longer over powers of a 2 (R) but over some different com- 
bination of powers of 1/R. 

2. Standard perturbative expansions 

2.1. Perturbative approach 

In this section we recall the standard perturbative method 
which is used to describe the early stages of gravitational 
dynamics in an expanding universe. For simplicity we fo- 
cus on a critical-density universe (i.e. f2 m = 1) but our 
results also apply to other cosmological parameters. Since 
we are interested in the large-scale behaviour of the den- 
sity field we adopt as usual an hydrodynamical descrip- 
tion. Thus, the matter distribution is approximated as a 
pressureless fluid with no vorticity. Besides, we only con- 
sider the linear growing mode. Then, one can look for a 
perturbative solution to the equations of motion for the 
density contrast #(x, t) and the divergence 0(x, t) = V.v of 
the peculiar velocity v(x, t) (indeed, within perturbation 
theory we can take the velocity field to be irrotational) . 
Here we introduced the comoving coordinate x. As shown 
in Goroff et al. (1986) the time dependence factorizes so 



that these expansions have the form: 

oo oc 

(5(x,i) = <5 (0 (x), 6>(x,t) = ^ho 1 w (x), (1) 



i=i 



i=i 



where a(t) oc t 2 ^ 3 is the scale-factor and a — da/dt is its 
time-derivative. The terms 5^ and 6^ do not depend on 
time and they are of order I over the growing linear den- 
sity field (5^^(x) = 5l(x). Thus, in the following we note 
<5l(x) the linear density field extrapolated up to z — (i.e. 
a = 1). In fact, since time plays no role (because of the 
factorization) we shall simply investigate the properties of 
the density contrast <5(x) today. The quasi-linear regime 
simply corresponds to the limit of small amplitude for the 
power-spectrum of the linear field <5l(x) (as described be- 
low this is usually different from the limit of large scale). 
In order to compute the perturbative expansions (|l|) it is 
convenient to work in Fourier space. Therefore, we define 
the Fourier transform of the density field as: 

/Hx 
^3 e ~' ik ' X ^ X )> ( 2 ) 

and similarly for all other fields. Then, substituting the 
exp ansions (|l|) into the equations of motion one obtains 
(see iGoroff ct al. (1986)|) : 



(qi, .., qz) with degree zero. They are symmetric over per- 
mutations of their arguments. These functions obey the re- 
cursion relations obtained in Goroff ct al. (1986) (see also 
Iain & Bertschinger (1994)). In addition, it is useful to 
note that the kernels F^ s ^(qi, .., qz) are proportional to k 2 



for I > 2, with k = qi + .. + qz, sec poroff ct al. (1986) 



This is a consequence of momentum conservation (e.g. 



Peebles (1974), Hogan et al. (1982)) 



Next, in order to derive the statistical properties of 
the density field <5(x) we must specify the properties of 
the linear density field $l(x) over which the perturbative 
expansions are built. In this article we investigate the case 
of Gaussian initial conditions and we take the linear grow- 
ing mode #l(x) to be a random Gaussian field. Therefore, 
its statistical properties are fully determined by its power- 
spectrum Pb{k) which we define by: 



(8 L {k)5 L {k')) = P L {k) 6 D (k + k') 



(5) 



Here we explicitly wrote the subscript "L" in order to 
recall that Pz,(&) is defined with respect to the field <5l(x): 
it is not the power-spectrum of the actual density contrast 
<5(x). We shall often consider power-law power-spectra of 
the form: 



P L (k)(xk n with -3<n<l, 



(6) 



which cover all cases of cosmological interest. Of course, 
the actual power-spectrum is not a pure power-law but it 
is usually sufficiently smooth (like the CDM model) to be 
approximated by a power-law over the range of interest. 

2.2. The non-linear power-spectrum 

From the standard perturbative approach recalled in 



Sect. [2.l| it is possible (at least in principle) to compute 
the various moments of the actual density contrast <5(x) 
up to any order I over the field <5l(x). In particular, 



some previous works ( 


e-g-, 


Jain & Bertschinger (1994), 


|3coccimarro & Frieman 


(1996) 


) have studied the non- 



linear power-spectrum P(k) of the density contrast i5(x), 
which is defined as in eq.(p|) by: 



(S(k)S(k')) = P(k) Mk + k'). 



(?) 



It is obtained by substituting the expansion (Q) into eq. ^ 
which yields up to order P 2 : 



«W(k) = ydq 1 ..dqzMk-q 1 -..-qz)F/ s) (q 1 ,..,qz) 
x *i(qi)..* L (q,) (3) 

and: 

{l) (k) = - J dqi..dqz Mk- qi -..-qz)G ; (s) ( qi ,..,qz) 

x S L ( qi ).Mm). (4) 

Here Sd is Dirac distribution while and are di- 
mensionless homogeneous functions of the wave vectors 



P(k) S D (k + k r ) = 




<5»(k) 




= (5«(k)5«(k')) + (<5( 3 )(k)(5 (1) (k')) 
+(S (2) (k)(5( 2 > (k')) + (S (1) (k),5 (3) (k')) + C{8 



(8) 



The various averages are computed from eq. and eq. (JsJ) 
us ing Wick's theorem. Such a stud y is detailed for instance 
in gcoccimarro fc Frieman (1996) where the expansion (|sj> 



is described as a sum over specific diagrams. Although 
these works were developed in Fourier space k it is conve- 
nient for our purposes to work in real space x, at least in 
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the intermediate steps. Thus, we can consider the expan- 
sion of the two-point correlation function £(x, x') defined 
by: 

e(xy) = (5(x)5(x')). (9) 
As in eq.(|^), we can introduce the expansion ([!]) through: 

£(x,xoh(x> (0 (x)) (f^'V)])- (io) 

In real space, the expression (||) for the terms 5^ (x) reads: 
5«(x)= J <ba..dxi F^frxx,..^) Sl^.Mxi) (11) 

(s) 

where the symmetric functions iy are the inverse Fourier 

transforms of the functions F$"\qx, ..,q/) introduced in 
eq.|): 

ifWi,--^) =|dk e i(k.x- q ,x 1 -..- q ,x I ) 



x F/ s) (qi, ..,q ; ) £ (k-qi - ..q ( ) 

.0), 



(12) 



The functions iy (x;xi, ..,x;) are invariant under global 
translations and rotations. Of course, since the initial con- 
ditions are homogeneous and isotropic the two-point cor- 
relation £(x, x') only depends on |x — x'|. The expansion 
( p^| ) is equivalent to the expansion (||). In particular, we 
have: 



ax) = / dk 



vk.x 



P(k), 



where we defined the function of one variable by 
£(xi,X2) = £(jxi — X2I). Then, using Wick's theorem the 
expansion ( |l 0[ ) can be written as a power-series over the 
two-point correlation A^(x, x') of the linear density field 
#£,(x), defined by: 



(6 L (x)5 L (x!)) = Al(x,: 



(14) 



Of course, the kernel A^(x, x') only depends on |x — x'| 
and the function Al(x) is the inverse Fourier transform 
of the linear power-spectrum Ph{k) introduced in eq.(||), 
as in eq.(O). As in eq.(p|), the expansion (poh reads up to 



(15) 



order A^: 

£(x,x') = (5«(x)5«(x')) + (5 (3) (x)<5 (1) (x')) 
+ {5 {2 \-x)5^ {■*!)) + (5«(x)5( 3 )(x')) +0(St) 

We give in Fig|l] a diagrammatic representation of the ex- 
pansion @ Such a geometric representation was already 
displayed in [Bernardeau (1992J 



Thus, the standard expansion (|l^) yields the two-point 
correlation £(x, x') as a perturbative expansion over the 
Gaussian random field 8^. This can be written as: 



£(x, x ') = ]T^W) 



(16) 



(1) 



(2) 



(3) 



(4) 



Fig. 1. Diagrammatic representation of the standard ex- 
pansion (|15| ) for the two-point correlation £(x, x') = 
(<5(x)<5(x')). The filled circles on the left represent the ac- 
tual density contrasts <5(x) and <5(x'). A white circle on 
the right with I points corresponds to <5")(x) and each dot 
represents a factor 8l in the expansion of 8^ . Each line 
stands for a factor A^ arising from the Gaussian average. 
We did not write multiplicity factors. 

where the term £W is of order A^. In particular, the 
lowest-order term (the "tree-level" contribution) is sim- 
ply: 

^)(x,x') = A L (x,x'). (17) 

It corresponds to the diagram labeled (1) in Fig.|. The 
next-to-leading term £( 2 ) is given by the diagrams (2) to 
(4) in Fig.[j], with appropriate multiplicity factors. Its ex- 
plicit expression is easily derived from eq.(lq). 



(13) 3. Small scale divergences 



The standard expansions built from the method re- 
called in Sect.|| have been used to obtain the cumu- 
lants of the density contrast at leading-order ( "tree- level" ) 



in 


Bcrnardeau (1992) 


and 


Bcrnardeau (1994) 


Then. 


Lokas et al. (1996) 


and 


Scoccimarro & Frieman (1996)| 



1=1 



rection). They noticed that for power-law linear power- 
spectra Pl with n > — 1 this subleading term actually 
diverges at large wavenumbers k — > 00. This ultravio- 
let divergence comes from small scales which are non- 
perturbative and which cannot be described by the hydro- 
dynamic equations of motion. This means that the stan- 
dard perturbative expansion breaks down for n > — 1. 
In fact, as we note below one can see that for n > — 3 
there always exists a finite order I beyond which some of 
the perturbative terms diverge. Therefore, the perturba- 
tive expansion is rather ill-defined. This is not surprising 
since for n > — 3 sufficiently small scales are far within 
the non-linear regime, which cannot be reached by per- 
turbative means. Indeed, perturbative approaches cannot 
go beyond shell-crossing (see for instance the discussion 
in jjaper I ). 

The high— k divergences arise from the integrations 
over the "small" branches contained within one of the 
two circles x or x' in Figjl]. This can be easily under- 
stood as follows. By looking at the two-point correlation 
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£(x, x') we select the large scale |x — x'| for the factors 
Al which arise from the vertical branches joining both 
circles. This provides a window function of scale |x — x'| 
which prevents any high— k divergence. By contrast, there 
is no such constraint for the factors which come from 
the small branches within one circle. Then, the integra- 
tion over the wavenumbers k associated with these small 
branches can lead to high— k divergences since they do 
not exhibit any window which would cut the contribu- 
tions from small scales^. This behaviour occurs for the 
diagrams (2) and (4) shown in Fig|l] for a steep power- 
spectrum n > — 1 but it only appears at higher orders for 
n < — 1. 

Thus, let us consider the contribution of diagrams of 
the form (2), with only one large branch of size ~ |x — 
x'| = R, one factor Sl inside the circle x' and 21 + 1 
factors Sl inside the circle x (Z = 1 in the diagram (2) 
shown in Fig.[l]). In Fourier space, such a diagram yields a 
contribution P\.2i+i{k) to the non-linear power-spectrum 
P(k) (k ~ 1/R) with: 

JW(fc) Mk + k') = (S L (k')S^ 2l+1 \k)) 

= /dqi..dq2/+i S D (k - qi - .. - q2i+i) 

x F^ +1 ( qi ,.. iq2l+1 ) (SlWlM.M^i+i))- (18) 
Using eq.(||) and Wick's theorem we obtain: 

Pl,2l+i(k) = (21 + 1)!! / dqi.-dq, P L {k)P L { qi )..P L {qi) 



(19) 



As we recalled in Sect.[2.1| the kernels F,' are propor- 



.0) 



tional to k 2 for Z > 2 (e.g., Goroff et al. (1986) ) because 
of momentum conservation. Therefore, introducing a high- 
wavenumber cutoff A (i.e. Pi(fc) = for k > A) we have: 



iW(fc) ~ P L (k) [ jj [A 3 P L (A)} 1 



(20) 



provided the integration over qi,..,q; diverges at large 
wavenumbers. For power-law power-spectra this condition 
corresponds to: 

(3 + n)l - 2 > 0. 



(21) 



Therefore, all such diagrams with I > 2/(3 + n) yield a 
contribution fc 2 Pi(fc) to the non-linear power-spectrum 
P(k). 

Here, we investigate the large-scale limit, which is dif- 
ferent from the early-time limit (if the power-spectrum 
is not a pure power-law). This is the limit of interest for 
practical purposes, where one tries to describe the density 
field on large scales which are still linear while the small 
scales have already entered the non-linear regime. Thus, 
we have shown above that all diagrams with only one large 
branch of size R and 21 + 1 factors Sl at point x or x' (and 

1 We would like to thank F. Bernardeau for discussions on 
these points. 



one factor Sl on the other side) in FigjT] yield a large-scale 
tail of order k 2 P L (k) for P(k), if I > 2/(3 + n). This corre- 
sponds to a term R~ 2 a 2 {R) for the non-linear two-point 
correlation £(i?) and variance £(-R). Here we introduced 
the rms linear density fluctuation a in a cell of radius R, 
volume V: 



a 2 (R) = (S 2 ) = 



f dxi dx 2 

/ — — A i(xi ,x 2 ), 



(22) 



where Sl.r is the linear density contrast within the spheri- 
cal cell of radius R. The non-linear variance £(i?) is simply 
given by: 



m = (4) 



(23) 



where Sr is the non-linear density contrast within the 
spherical cell. Therefore, we see that the large-scale non- 
linear two-point correlation £(i?) cannot be expanded as 
a series over powers of the linear variance a 2q , with q> 1. 
Indeed, beyond some order q there appears a term of or- 
der R~ 2 a 2 (R) arising from the diagrams described above. 
Note that diagrams with one branch of size R and with I 
(resp. I') factors Sl at the point x (resp. x') with both 
I > 1 and /' > 1 yield a contribution smaller than 
R~ 2 a 2 (R) by additional powers of a 2 or R~ 2 (or an entire 
power of a 2 ). On the other hand, the contribution from 
diagrams with I branches of size R (e.g., I — 2 for the dia- 
gram (3) in Fig.|l|) exhibits at least a factor a 21 . Therefore, 
the first "non-regular" contribution to the non-linear two- 
point correlation (i.e. not of the form a 2q ) is the term 
R~ 2 a 2 {R) described above. 

The order q + of the last "regular term" is given by the 
constraint q+(n + 3) < 2 + (n + 3) which yields: 



(24) 



where [..] is the entire part. Therefore, we see that for all 
power-spectra the expansion of the non-linear correlation 
£(i?) as a series over powers of a 2 breaks down beyond 
a finite order q + . This order q + is larger for smaller n 
which corresponds to initial conditions with less small- 
scale power. In particular, for n > — 1 the only regu- 
lar term in the expansion of £ (or P(k)) is the leading 
order-term a 2 since q + = 1. This means that the so- 
called "one-loop term" shows a small-scale divergence. 
This agrees with the results of the previous studies (e.g., 
[Lokas et al. (1996) , Scoccimarro & Frieman (1996)). For 
n < — 1 the break up of the expansion over powers of a 2 
only occurs at a higher order. Thus, in this respect the in- 
dex n = — 1 does not correspond to a significant change of 
the physics. Here, we note that two different cases may oc- 
cur with respect to the new contribution R~ 2 a 2 (R) which 
arises from an infinite series of diagrams. 

Firstly, the high-wavenumber cutoff A (in case it is set 
by the initial shape of the linear power-spectrum Pj,(fc)) 
may still be within the linear regime. Then, the pertur- 
bative expansion makes sense and one may still be able 
to write a perturbative expression for the large-scale non- 
linear two-point correlation £(R). However, in addition to 
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terms of the form a 2q one must take into account the 
new contributions like R~ 2 a 2 (R) which arise from par- 
ticular classes of diagrams. This requires some partial re- 
summations over infinite series of diagrams which may 
still lead to finite results (a definite statement requires 
an explicit calculation). If we also expand over powers of 
A 3 (A) ~ cr 2 (l/A) Cl we obtain a finite number of di- 
agrams, up to the required orders a(R) 2qi and a(l/A) 2q2 . 
Note that in any case this implies that the physics at 
scale R also depends on the shape of the linear power- 
spectrum at the much smaller scale 1/A. In fact, this is 
not totally surprising. Indeed, although we could expect 
the large-scale behaviour to be "largely" independent of 
the small-scale physics, this decoupling is not exact. That 
is, by choosing different cutoffs A we clearly obtain dif- 
ferent density fields and this must also show up, to some 
degree, in the exact properties of the density field at large 
scale R, even if R ^> I/A. Thus, our analysis shows that in 
the limit R <C 1/A <C Rq, where Rq is the comoving scale 
associated with the transition to the non-linear regime 
(defined for instance by cr(i?o) = 1) ; the dependence on 
A does not show up in the leading-order term of £ (which 
is merely a 2 ) but it can be traced in higher-order correc- 
tions. This clearly agrees with intuitive expectations. 

Secondly, the cutoff A may already be within the non- 
linear regime. This is actually the case of cosmological 
interest where the smaller scales have usually turned non- 
linear at the time of interest. Then, it is clear that the 
perturbative approach fails and rigorous results require a 
non-perturbative method. Thus, the perturbative contri- 
butions which are dominated by the cutoff A, as in eq. (p(j|) , 
cannot yield correct results since large wavenumbers in the 
range 1/Rq < k < A cannot be described by perturbative 
means. Indeed, on these small scales one must take into 
account shell-crossing, which is out of reach of any ap- 
proach which is based on the hydrodynamical description. 
One might still be able to estimate the order of magni- 
tude of the terms like eq.(pO|) by setting A = l/Ro, in the 
sense that the non-linear dynamics itself may generate the 
required cutoff to make the diagrams finite. However, it 
is clear that this heuristic procedure requires some non- 
perturbative calculations in order to get any convincing 
support. Besides, the divergence of subleading terms casts 
some doubts on the validity of the perturbative results ob- 
tained at lower orders. We shall come back to these points 
in Sect.[| and Sect J5| below. 

Finally, we note from eq.(pO|) that the degree of di- 
vergence of the terms encountered in the perturbative ex- 
pansion can be made arbitrarily large since it increases 
at higher orders, as A^ 3+n ^ -2 . Therefore, this is a non- 
renormalizable theory and standard renormalization pro- 
cedures cannot be applied (e.g., Zinn- Justin (1989)). 



4. Steepest-descent method 



In a previous article ( paper II ) we developed a steepest- 
descent method to study the properties of the density field 
in the quasi-linear regime. Here we describe how the anal- 



ysis performed in the previous section from the standard 
perturbative approach translates within this new frame- 
work. We study the statistics of the density contrast 5 R 
within a spherical cell of comoving radius R, volume V: 



5(x). 



(25) 



To this order, we introduce the generating function ip(y) 
defined by the relation: 



■4>{y) = ( e -a<W<7 2 (A)) 



dS R e- ySR/a2{R) V{S R ). (26) 



The average (..) in eq.(|26j) denotes the mean over the ini- 
tial conditions, which are defined by the linear growing 
mode (Sl(x). The last equality defines the probability dis- 
tribution function (pdf) Vidji) of the density contrast. 
Thus, ifi(y) is also the "rescaled" Laplace transform of 
"P((5r). It also yields the moments of the density contrast 
5 R through the expansion: 



{-y) q (s q R ) 



T 2 'l 



(27) 



9=0 



It is convenient to introduce also the generating function 
(p{y) of the cumulants defined by: 



My) = e-*WI°\ 

which exhibits the expansion: 



f(y) 



9=2 H 



i-y) q (S R )c 



,2(9-1) 



(28) 



(29) 



since we have (Sr) c — 0. Thus, the term of order y 2 of the 
expansion over y of the generating functions ip(y) or ip(y) 
yields the non- linear variance £ = (S R ) C - 

Then, as shown in paper II the generating function 
ip(y) can be expressed as a path-integral through the re- 
lation: 



</%) = (DetA^ 1 



-l\l/2 



[<Wi(x)] e 



-S[6 L ]/v 2 (R.) 



where we introduced the action S[Sl]'- 



S[S L ] = y 5 R [5 L 



a 2 (R) 



(30) 



(31) 



Here <5r[5l] is the non-linear density contrast 5 R within 
the spherical cell V obtained for a linear density field 
<5l(x). The relation ( |30| ) is merely the explicit expression of 
the Gaussian average over the linear growing mode <5l(x) 
which enters the definition of the generating function ip{y) 
in eq. (p6|) . 

The action S[5l] is independent of the normalization 
of the power-spectrum Ph{k) since A/, oc a 2 . Then, it 
is clear that the path- integral in eq.(|30|) is dominated 
by the minimum of the action S in the limit a — ► 
for a fixed y. Indeed, the contributions from other 
states Sl are exponentially damped relative to this point. 
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Moreover, the steepest-descent approximation becomes 
exact in this limit. Here we must point out that this actu- 
ally corresponds to the "early-time" limit and not to the 
"large-scale" limit. Indeed, the shape of the linear power- 
spectrum Ph{k) (or of the kernel A/,) is held fixed while 
its amplitude goes to zero. 

Next, let us note S y the saddle-point of the action S, 
for a given y, and S y = S[S y ] the value of the action at 
this point. Expanding the action around the point S y we 
write: 



S[Sy + S' L 



Sy 



12 



(32) 



where we introduced the symmetric kernels W q defined by 
the functional derivatives at the point 5 V : 



Wq(xi, ..,Xg) = 



5^S 



«(«L(xi))..<y(5£(x,)) 

and we used the short-hand notation: 

W q .5'l = J dxi..dx,W,(x 1 ,..,x fl )^(xi)..^(x 9 ). 



(33) 



(34) 



in Sectj^. For instance, as shown in paper II, for power- 
spectra with n < this method yields the correct slow 
decline at large densities of the pdf V{8b) (shallower than 
a mere exponential) which cannot be rigorously derived 
from the usual perturbative approach. This corresponds 
to the case where the radius of convergence of the se- 
ries expansion ( p9| ) is zero. The reason why it is possi- 
ble to go beyond the standard perturbative results is that 
the saddle-point S y is obtained in a non-perturbative way 
(this is made possible by the fact that the action S[Sl] 
is spherically symmetric and we know the explicit solu- 
tion of t he spherical dynamics b efore shel l-crossing, see 
Ipapcr II ). Besides, as shown in |paper IV the steepest- 
descent method can also be applied to the non-linear 
regime where it provides the tails of the pdf V(8r). 

For practical purposes, it is convenient to write eq.([37|) 

as: 



As noticed in |paper II , for small y the kernel W% is definite 



<P(V) =S y + — Trln(l + A y ) + .. 
where we defined the matrix A y by: 



A B (xi,x 2 )=^ y*dx A L (xi,x) 



$($ £ (x))*(fo(xa)) 



positive since W%{y — 0) = a A^ , see eq.(|31|). Then, we 
can use the steepest-descent method in the limit a — > 
0. The quadratic term (W 2 .£f )/a 2 in the exponential in 
eq.(30) provides a cutoff of order a for S' L so that the next 
term (W^.S'l) / a 2 is of order a — ► 0. This confirms that 
at leading order the generating function ip(y) is given by 
the Gaussian integration around the saddle-point S y . This 
yields: 



The value S y of th e action at the saddle-point 5- 



(38) 



(39) 



was 



derived in |paper II . It yields the leading-order contribu- 
tion to the generating function (p(y) and to the cumulants 
(SV) C , through eq.(29). In this article we are interested 
in the subleading corrections. In particular, the next-to- 
leading contribution is given by the term of order a 2 writ- 
ten in eq. ( |3e] ) . In order to derive the correction associated 
with a cumulant of order q we simply need to expand the 
_ , 2 logarithm ln(l + A y ) as a power series over y, up to order 

ip(y) = (DetA^ 1 ) (DetAf a )~ 1/2 e ~ s «/ CT ^ + ... (35) y i. Here, we shall only consider the cumulant of order 2 



where the matrix M y is the Hessian of the exponent at 
the point 5 y : 



(that is the power-spectrum or the two-point correlation) 
so that we can write eq.(38) as: 



Mj,(xi,x 2 ) 



S 2 (S/a 2 



y s 2 (S R ) 



y 2 (S 2 R )c 



V 2 v 2 , 
-?-+— TrA„ 



-TrAl 
2 y 



.(40) 



a 2 (R) 5(* L (xi))*($£(x 2 )) 



where we used S y = 



-y 2 /2 + <D(y 3 ) (see paper II) and the 



a -1/ \ /op\ fact that A y is of order y (see eq.(|39|)). Note that we must 



Higher-order corrections to ip(y) (i- e - higher powers of a 2 ) 
are obtained by expanding the terms of order q > 3 of the 
action S in the exponential in eq. (|30|) . In fact, the expres- 
sion (|35|) already gives the next-to-leading order term of 
the cumulants (<5fj) c . Indeed, from eq.(p8|) we have: 



check that we recover (8r) c = 0. We shall also recover the 
results of Sect.||. 

Now, we need the expression of the matrix A y up to 
order y 2 . Therefore, we must obtain the second-derivative 
which appears in eq.(p9|) up to order y. Moreover, as seen 



in paper II the saddle-point S v is of order y in the limit 



cp(y) = -cr 2 lnV> = S y + — lnDet(Ai.M y 



(37) 



This provides the generating function ip(y) as an expan- 
sion over a 2 . As explained above, the higher-order terms 
in eq.(|37]) come from the expansion of the exponential of 
the action beyond its quadratic term in eq.(|3p|). Here we 
must point out that although the steepest-descent method 
yields an expansion over a 2 it can actually go beyond 
the reach of the standard perturbative approach described 



y — * 0, see also eq. (J45|) below, so that we can use the 
expansion ( [fl| ) to write up to order y: 

S(x)[S v +6' L ] = J d Xl dx 2 F 2 (s) (x; Xl ,x 2 ) &(xi)^(x a ) 
+3 J dx 1 ..dx 3 F 3 (s) (x;x 1 ,..,X3)^(x 1 )^(x 2 ) ( 5,(x3)(41) 

where we only kept the terms of order S' 2 . Hence we have 
up to order y 2 : 

4,(xi,x 2 ) = J / y/ dx ' A L (x 1 ,x')F 2 (s) (x;x',x 2 ) 
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~j Y J dx'dx" A i (x 1 ,x')^ (s) (x;x',x",x 2 ) ( 5 3/ (x") 
,. V (42) 



From eq.(40) we obtain for the cumulant of order 1, 
up to order a 2 : 

(Sr) c = J v yJ dx'dx"A i (x",x')F 2 (s) (x;x',x"). (43) 



Using the standard result (e.g., [Peebles (1980)[ ) 



f'W(k 1 ,k a ) = ~ + ~ 



5 1 ki.k 2 1 ki.k 2 2 /ki.k 2 



7 2 kf 2 k 2 7 V hk 



(44) 



it is easy to show that we recover (<5r) c — 0. Then, from 
paper II we can write the saddle-point 8 y (x) up to order 
y as: 

<5,(x) = -4/ ^ A L (x,x')+0(2/ 2 ). (45) 

This yields for the cumulant of order 2, up to order a 4 : 

(S 2 n)c = a 2 

+ 6 J ~^T~ J dx l'- dx 3 •F 1 3 (S) ( X ; X 1: X 2,X 3 ) 

x Al(xi,x 2 )Al(x 3 ,x') 

+ 2 y dxd^ y dxi .. dX4F ( S ) (x;XijX2)F M (x /. X3;X4) 



x A L (xi,x 3 )A L (x 2 ,x 4 ) 



(46) 



Thus, we recover the result we would obtain from the stan- 
dard expansion displayed in Fig.[lJ This is not surprising 
since we used the same expansion ( |IT|) to write the fluc- 
tuations of the functional <5r[<5l] in eq.(|4^). 

Of course, as in Sect.||, in the case of power-law power- 
spectra we must add a high-wavenumber cutoff A in order 
to obtain finite results. Note that within this steepest- 
descent method we consider the "early-time" limit where 
the amplitude of the power-spectrum goes to zero, and 
not the "large-scale" limit. Therefore, <r 2 (l/A) is of the 
same order as cr 2 (R) (they are both proportional to the 
amplitude of the linear power-spectrum, whose shape is 
held fixed) and the expansion obtained from eq. (|3(]) is over 
powers of both cr 2 (l/A) and a 2 (R), which are treated on 
the same footing. In the limit A — > oo we again encounter 
divergent quantities, as in Sect|| 

In order to go beyond these divergences we should use 
a non-perturbative expression for the functional <$r [#£,]. 
In particular, perturbative expressions like eq.(|4l|) for the 
fluctuations of the non-linear density contrast Sr around 
the saddle-point 8 y are not sufficient. However, we note 
that the saddle-point 8 y and the value S y of the action 
at this point are obtaine d in a n on-perturbative man- 
ner. Indeed, as shown in paper II this calculation does 
not rely on a perturbative expansion for the density field 
(5(x). Therefore, the quantity S y is the exact minimum of 
the action even if we use the exact non-linear ex- 

pression for the functional 8r[6l\- Hence the leading-order 



term S y for the generating function ip{y) is correct and it is 
not affected by the divergence of subleading terms beyond 
some finite order. From the expansion ( p9| ) this also means 
that the value obtained for the cumulants at leading-order 
("tree-level") is exact. 

This agrees with the fact that numerical simulations 
are seen to match the predictions of leading-order per- 
turbation theory which is consistent with eq.( |3^) for the 
see 



cumulants (8 q 



H/ci 

This result could also be expected within the frame- 
work of the usual perturbation theory. Indeed, as shown 
in Bernardeau (1992,1994), at leading-order in the limit 
a — ► the cumulants only involve "tree diagrams" with 
no "small internal branches" which were at the origin of 
the small-scale divergences which appeared for instance in 



Bernardeau (1994) and paper II 



As discussed below in Sect. 5.2, we can actually expect 
the results of perturbation theory to be valid up to the 
order q + obtained in eq.(|24|), that is until they are finite. 

5. Comparison with two simple dynamics 

We have seen in the previous sections that for power-law 
power-spectra Pb{k) ex k n one always encounters some 
divergent quantities beyond some finite order in a per- 
turbative expansion over a 2 . However, in the quasi- linear 
regime the results from numerical simulations agree very 
well with the leading-order terms obtained by such an- 



alytical means (e.g., |Bernardeau (1994)[ ). Of course, nu 



merical simulations always have a small-scale cutoff due 
to finite resolution effects but it has been checked that 
the behaviour they obtain converges and becomes inde- 
pendent of this cutoff (i.e., the limit A — > oo appears to 
be well-defined). This suggests that the divergence of these 
subleading terms does not spoil the predictions obtained 
from the calculations of the lower-order terms. 

One may then argue that the non-linear dynam- 
ics itself could provide a small-scale cutoff Rq. Indeed, 
at scales smaller than Rq which marks the transi- 
tion to the non-linear regime (a(Ro) — 1) perturba- 
tive approaches break down so that the integrals ob- 
tained in the previous sections do not make sense at 
large wavenumbers k > 1/Rq. Therefore, one should 
cut the integrations at A = l/i?o and the contribu- 
tions of higher wavenumbers should be derived from a 
fully non-perturbative method. Then, this could pro- 



vide finite results (see also 


Jain & Bertschinger (1996) 


|3coccimarro & Frieman (1996) 


. for discussions on this 



point). Unfortunately, a more detailed statement requires 
explicit non-perturbative calculations which remain to be 
done. In order to get some insight into this problem, it 
could be useful to study a somewhat simpler dynamics as 
described below. 

5.1. The Zel'dovich approximation 

For instance, we may try to investigate these small-scale 
divergences within the framework of the Zel'dovich ap- 
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proximation ( Zel'dovich (1970)). Unfortunately, we shall 
see below that these divergences actually disappear in this 
simpler dynamics, if it is to make any sense. Indeed, as 
shown in Grinstcin fc Wise ( 1987)[ , within the framework 
of the Zel'dovich approximation (ZA) the symmetric ker- 
nels which enter the expansion (|IJ) through eq. (||) are 
given by: 



(ZA): F['\q 1 ,..,q i ) = ^ 



Therefore, we have F^^Qs., qi, — qi, .., q;, — q;) oc k' M and 



1 k.qi k.q/ 



Qi Qi 



(47) 



eq.(20) becomes: 



(ZA): P ll3I+1 (fc) ~ P L (k) (£) [A 3 P L (A)]' 
oc A(" +1 )' 



(48) 



provided the integral is dominated by the cutoff A. We 
see that this constraint requires n > — 1. However, for 
such power-spectra we cannot use the Zel'dovich approx- 
imation because the linear peculiar velocity field v/, (x) is 
not well-defined. That is, high wavenumbers k provide a 
divergent contribution to the unsmoothed velocity (but of 
course the mean velocity over some finite spherical cell V is 
finite). Therefore, we cannot study this dynamics without 
imposing some high-wavenumber cutoff A. Note that one 
sometimes uses the non-linear scale i?o as a small-scale 
cutoff, which is refered to as the "truncated Zel'dovich 



approximation" (e.g., Coles et ah (1993) ). On the other 
hand, for n < — 1 there are no small-scale divergences 
at all: all diagrams are finite. This allows one to write a 
perturbative expansion for the non-li near two-point corre - 
lation £(i?) where all terms are finite ( Bharadwaj (1996) ). 
Therefore, the Zel'dovich and the exact gravitational dy- 
namics exhibit rather different behaviours in this respect. 
Thus, we cannot use the simple Zel'dovich approximation 
to investigate the small-scale divergences. However, we can 
use the adhesion model which still makes sense for n > — 1, 
as described below. 



5.2. The one-dimensional Burgers equation 

The adhesion model allows one to mimic the gravi- 
tational sticking of particles after shell-crossing which 
is not described by the Zel'dovich approximation, see 
Gurbatov et al. (1989)| . This provides a reasonable de- 



scription of the large-scale skeleton of gravitational struc- 
tures (e.g., Vergassola et al. (1994)). To do so, one sim- 
ply adds a vanishing viscosity to the Zel'dovich equations 
of motion. Then, as shown in Gurbatov et al. (1989) (see 



also Vergassola et al. (1994)) the evolution of the velocity 
field (after suitable rescaling) is described by the stan- 
dard Burgers equation. The infinitesimal viscosity pre- 
vents shell-crossing as particles stick together after colli- 
sion. Then, contrary t o th e simple Zel'dovich approxima- 
tion discussed in Sect.ST we can study initial conditions 
with n > — 1. Indeed, the non- linear dynamics itself makes 



the velocity field finite at any time t > 0. Here we are only 
interested in the adhesion model for illustrative purposes 
hence we consider the standard one-dimensional Burgers 
equation: 

^ (49) 



dv 



dv 

dt dx dx 2 

in the inviscid limit v — > 0. Moreover, we investigate the 
case of a white-noise power-spectrum P V(S (k) for the initial 
Gaussian velocity field vq(x) at t = (i.e. P Vo (k) is flat): 



(v {k 1 )v (k 2 )} 



Pv (ki) S D (k 1 +k 2 ) 
r S D (ki + k 2 ). 



(50) 



Here the length xq and the time tg arc parameters which 
merely define the normalization of the power-spectrum 
P vo (k). Then, the two-point correlation of the initial ve- 
locity field is: 



A Vo (xi,x 2 ) = (vo(xi)v (x2)) 



2;r .r„ ( ^ ) ,) v (.r l - - r 2 I. 



(51) 



The kernel A Vo fully defines the statistics of the Gaussian 
random field Vq(x). Then, the velocity field v(x,t) at time 
t > is obtained by solving Burgers eq.([49"|) with the ini- 
tial condition Vq(x). This classical problem is studied in 
details in Burgers (1974)| (se e also Gurbatov et al. (1989) 
and Vergassola et al. (1994)). In particular, the exact so- 
lution v(x, t) can be obtained through the standard Hopf- 
Cole transformation. 

Starting at t — with a uniform density distribution 
it is easy to see that at any regular point x the density 
contrast 5{x,t) is given by: 

S(x) = -ig. (52) 

Therefore, in order to obtain the statistics of the den- 
sity field we merely need the properties of the velocity 
field. First, we can investigate a perturbative approach 
in order to derive v(x,t), in a fashion similar to the pro- 
cedure detailed in Sect.|2| for the gravitational dynamics. 
From eq. ( |49"|) we obtain at linear order in the inviscid limit 
(y — ► 0) the simple relation dv/dt = 0. Thus, the linear 
solutions Vl(x, t) and 8l{x, t) are: 

,dv 



v L (x,t) = v (x), 8 L (x,t) = -t 



dx 



(53) 



Then, as in eq. (|l]) we can look for a perturbative solution 
v(x,t) and 5{x, t), written as an expansion over the ini- 
tial velocity field vq{x) (i.e. over the linear growing mode 
5l(x, t)). Shell-crossing does not appear in the perturba- 
tive approach (e.g., Bharadwaj (1996), paper ]) so that in 
the limit v — > we actually recover the Zel'dovich dynam- 
ics (i.e. the r.h.s. in eq.p9|) disappears as v -— > within 
perturbation theory). Then eq.(]49]) reads dv/dt = and 
the distribution function f(x,v;t) is simply: 

f(x, v; t) — po 8o(v — vq(s)) with x(s,t)=s + tv, (54) 
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where po is the initial uniform density. Next, expanding 
the Dirac in cq.( 54 ) and using p(x,t) — J f{x,v;t)dv we 
obtain as in Grinstcin & Wise (1987): 



5(x,t)=J2 



(-!)< 



1=1 



Ql 

dx 1 



[v (xY] . 



(55) 



In Fourier space this yields: 
°° 1 f 

s ( k ) = X] 7? / d 9i-- d « S o(k -qi- ..- qi) 



i=i 



qi--qi 



SL(qi)--S L (qi)- 



(56) 



Hence we recover the one-dimensional form of eq.(K 



for the kernels F; (<7i, .., qi). From eq.(|53|) the power- ((<5^ 
spectrum P^k, t) of the linear density field Sl(x, t) is sim- 
ply: 



P L (k,t) = t 2 k 2 P Vo (k) oc k 2 
so that eq.(|4"8|) becomes: 



21 



Pi.2i+i(k) ~ P L {k) [ j ) [AP L (A)} 1 oc A 1 



(57) 



(58) 



Thus, we see that we again encounter divergent quantities 
in perturbation theory, as in Sect.|| for the gravitational 
dynamics. Note that these divergences already appear at 
the next-to-leading order, for I = 1 in eq. (p8|) . 

In fact, for a white- noise initial velocity the density 
contrast within a top-hat of radius R is not well-defined, 
even at linear order. Indeed, the initial velocity field vq(x) 
is not a smooth function but a distribution. Thus, we 
study the mean density contrast over a Gaussian win- 
dow of size R: 



Ax 



2nR 



S(x). 



(59) 



The Gaussian cutoff ~ e~( kR ^ I 2 at high wavenumbers en- 
sures that ((Sr) 2 ) is well-defined at linear order. By con- 
trast, a top- hat window exhibits an extended power-law 
tail at large k. Thus, in order to really probe the physics 
at scale R it is better to use the Gaussian window. Using 
eq.(52) we can write: 



and: 

<(^) 2 > 



l^R? 



dx e -* 2 /(2* 2 ) x 



dx 



-x 2 /(4R 2 ) 



1 



X 2 \ 

2i? 2 ) 



(60) 



A„(a;).(61) 



Here we introduced the two-point velocity correlation 
A„(x) = (v(r)v(r + x)}. At linear order we have v — vl = 
vq, hence A v (x) = A v0 (x) and using eq.(51) we obtain: 



; / t 

T 



(62) 



Note that small scales turn non- linear first, following a hi- 
erarchical process as for the usual cosmological scenario of 
structure formation. Then, a perturbative approach would 
give the non- linear variance ((S^) 2 ) as an expansion over 



powers of <7q. However, as noticed from eq.(58) the per 



turbative method already fails at the next-to-leading level. 

The advantage of the Burgers equation (49) is that it 
can be solved explicitly through the Hopf-Cole transfor- 
mation. This allows one to derive numerous exact results 
in a non-perturbative way (e.g., Burgers (1974)| and ref- 
erences in Vergassola et aL (1994)). In particular, we can 
obtain the variance ({S'^) z ) as an expansion over 1 / R in 
the large-scale limit R — > oo by expanding the exponential 
in eq.(pl[). At leading order we get: 



t 2 



40Fi? 3 J- 



dx A„(a;). 



(63) 



Next, we note that the integral J = J dxA v (x) is ac- 
tually independent of time (i.e. it is an invariant, see 
p3urgcrs (1974) ), as can be seen from eq.([49|). Using 
eq.(|5l|) we obtain: 



J — J dx A v (x) — J dx A v q(x) = 2tt 



Substituting this result into eq.(p3[) yields: 

3 / t n - 
to 



(64) 



(65) 



Thus, at leading order in the large-scale limit R — > oo we 
recover the result (62) of perturbation theory. In a similar 
fashion, the next-to-leading order term is: 



<(^) 2 > (2) 



-3t 2 



160Fi? 5 J_ 



dx x 2 A v (x) oc R 



(66) 



Indeed, 
ishes as 



the 



two-point correlation 

-\x\ 3 



A v (x) van- 
at larg e scales |x| — * oo, see 
|Frachebourg fc Martin (2000) , so that the integral 



converges. Thus, we can see that the term ((Sr) 2 )^ 
decreases more slowly than <Tq at large scales. This 
explains why the perturbative approach must break down 
at the next-to-leading order. In fact, we see that the 
expansion of ((S^) 2 ) is of the form: 



((<*: 



G\2\ 



3-2(;-i) 



(67) 



i=i 



This expansion is not over powers of Oq. In fact, beyond 
the order q + = 1 where perturbation theory is valid there 
are no more traces of the naive perturbative expansion 
since no powers of <Jq appear. On the other hand, we note 
that the first "non-regular" term (i.e. which is not a power 
of Gq) of the exact expansion (p7|), obtained in cq.(pq), is 
actually of the form (|5S|) , with 1 = 1. Thus, the non-linear 
dynamics has introduced a cutoff A = 1/Rq, where Rq 
is the scale which marks the transition to the non-linear 
regime. This makes the first divergent term ( |58|) finite. 



10 



P. Valageas: Dynamics of gravitational clustering V. Subleading corrections in the quasi-linear regime. 



However, this is not the whole story since the form of the 
exact large-scale expansion ( |67| ) beyond this order cannot 
be infered from the naive perturbative theory (|p) and 
(|56|), which would also give terms of the form a G 1 R~ 21 " 2 
which do not appear in the exact result (|67|). 

We can reasonably expect the previous results ob- 
tained for the Burgers equation to apply to the case 
of the gravitational dynamics described in Sect.^J, al- 
though a definite statement would require an exact 
calculation. Indeed, the Burgers dynamics has been 
seen to provide a good descri ption of gravitationa l 



clustering on large-sca les (e.g., |Gurbatov et al. (1989) 



Vcrgassola ct al. (1994)| ). Moreover, the physics involved 
in the divergences of the perturbative expansions appears 
to be similar. Thus, in both cases the break down of 
the perturbative approach is related to its failure to take 
into account shell-crossing. Moreover, when particles cross 
each other caustics appear where the density is infinite and 
it is clear that perturbative methods cannot go beyond 
these points. Then, in the adhesion model particles stick 
together while for the exact gravitational dynamics parti- 
cles actually cross each other and build "virialized" objects 
of finite thickness stabilized by a non-zero velocity disper- 
sion. Nevertheless, in both cases the dynamics provides 
a cutoff at the non-linear scale Rq. For instance, while 
the typical displacement of particles diverges in perturba- 
tion theory, because of the contribution of high wavenum- 
bers, it is clear that the actual contribution of small scales 
(R < Ro) is finite and of order Rq. 

Thus, by analogy with eq.(|6^) we can expect the 
perturbative results obtained through the standard ex- 
pansions ([!]) to be valid in the large-scale limit up 
to the order q + derived in eq.(|24|), that is until they 
are finite. This agrees both with the analysis described 
in Sect.fl which shows that the leading-order terms 
must be correct and with the calculations performed in 



Scoccimarro fc Frieman (1996) which find that numerical 



simulations match the predictions of perturbation theory 
up to the next-to- leading order for n < —1, when they 
are finite. Moreover, by analogy with eq.(|5e]) and eq. (|66|) , 
we can expect from eq.(^0]) that the first "non-regular" 
term (i.e. not of the form a 2q ) in the expansion of £(i?) 
over powers of l/R in the large-scale limit is of the form 
R~ 2 (T 2 (R). Finally, we note that in order to get meaning- 
ful results it may be useful to consider a Gaussian window 
rather than a top-hat which exhibits an extended tail at 
high wavenumbers. 

6. Conclusion 

Thus, in this article we have investigated the properties of 
the standard perturbative expansions which arise from the 
dynamics of gravitational clustering. First, we have shown 
that for hierarchical scenarios with no small-scale cutoff 
perturbative theory always breaks down beyond a finite 
order q + as divergent terms appear. Moreover, the degree 
of divergence increases with the order of the perturbative 
terms so that standard renormalization procedures do not 



apply. At first sight, the divergence of these subleading 
terms may cast some doubts on the predictions of pertur- 
bative theory obtained at lower orders. However, we have 
shown that the leading-order results can also be derived 
through a steepest-descent method which does not make 
use of these perturbative expansions. Therefore, the pre- 
dictions obtained at leading-order for the cumulants of the 
density contrast are correct. This agrees with the results 
of numerical simulations which match these calculations 



(e.g., Bernardeau (1994)). Then, we have shown that the 
simpler one-dimensional Burgers equation exhibits a sim- 
ilar behaviour. This analogy suggests that the results of 
perturbation theory are valid up to the order q + (i.e. until 
they are finite) and that beyond this order the large-scale 
expansion of the two-point correlation is no longer over 
powers of cr 2 (R) but over some different combination of 
powers of l/R. 

In practice the linear power-spectrum shows a cut- 
off at a very small scale R c , like for the CDM model 



(Peebles (1982)). This ensures that all integrals encoun- 
tered in perturbative expansions are finite. However, it is 
clear that the results of perturbative theory break down 
when the contribution of non-linear scales becomes non- 
negligible. Thus, our analysis also applies to the cases 
where R c < Rq, where Rq is the scale which marks 
the transition to the non-linear regime. Then, to go be- 
yond the order q + one needs again to devise new non- 
perturbative methods. On the other hand, if R c > Ro the 
standard perturbative expansions are sufficient, but this 
is not the case of cosmological interest. We can also note 
that although we considered a critical density universe so 
far, our analysis also extends to arbitrary cosmological 
parameters. 

Finally, we can point out that the leading-order terms 
are quite sufficient to describe the quasi-linear regime, up 
to a ^ 1, as seen in Bernardeau (1994) and [paper II . 
Therefore, the break-up of perturbative expansions be- 
yond a finite order is not a serious problem. Moreover, as 
discussed in paper 1 (see also Bharadwaj (1996)) pertur- 
bative expansions do not capture any effect arising from 
multi-streaming. Hence, in order to derive convincing re- 
sults in the early non-linear stages a <; 1 where a signifi- 
cant fraction of matter has experienced some shell-crossing 
(let us recall that a spherical density fluctuation collapses 
to a singularity as soon as 8l ~ 1.68) one should in any 
case use some non-perturbative approach. This still re- 
mains an open problem. However, we can note that some 



results in this direction have been obtained in paper IV 



for the rare-event tails of the probability distribution of 
the density contrast. 
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